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We consider a topologically non-trivial flat band structure in one spatial dimension in the presence 
of nearest and next nearest neighbor Hubbard interaction. The non-interacting band structure is 
characterized by a symmetry protected topologically quantized Berry phase. At certain fractional 
fillings, a gapped phase with a filling-dependent ground state degeneracy, and fractionally charged 
quasi-particles emerges. At filling 1/3, the ground states carry a fractional Berry phase in the 
momentum basis. These features at first glance suggest a certain analogy to the fractional quantum 
Hall scenario in two dimensions. We solve the interacting model analytically in the physically 
relevant limit of a large band gap in the underlying band structure, the analog of a lowest Landau 
level projection. Our solution affords a simple physical understanding of the properties of the gapped 
interacting phase. We pinpoint crucial differences to the fractional quantum Hall case by studying 
the Berry phase and the entanglement entropy associated with the degenerate ground states. Finally, 
the symmetry protected nature of the interacting phase is demonstrated by explicitly constructing 
a gapped interpolation to a state with a trivial Berry phase. 

PACS numbers: 03.65.Vf, 73.43.-f, 72.15.Nj 






I 

o 
o 



> 

\o 
en 

o 
en 

> 



I. INTRODUCTION 

Interacting topological phases that can not be un- 
derstood at the level of non-interacting models have 
attracted continued interest since the discovery of the 
fractional quantum Hall (FQH) effeclPt^l. pQH states 
can be observed in partially filled Landau levels, i.e., in 
systems that would be metallic in the absence of correla- 
tions. At certain fractional fillings, the huge phase space 
for interactions provided by the macroscopic degeneracy 
at the Fermi level then co nspir es with the non-trivial 
topology of the Landau levepl^ to yield a gapped state 
which is topologically distinct from any non-interacting 
two-dimensional (2D) insulator. To make this distinction 
more precise, the notion of topological order has been 
introduced by WerP. The simplest FQH state can be 
observed at i/ — ^ filling. From a phenomenological 
viewpoint its key differences from the integer quantum 
Hall (IQH) state"* ^ ^ observed in a completely filled 
Landau level are the following: (i) The Hall conduc- 
tance axy representing the topological invariant of the 
IQH state assumes a fractional value, more concretely 
a^y = \\- (ii) If periodic boundary conditions are 
imposed, the system exhibits a threefold degenerate 
ground state for v = ^. (iii) The elementary excitations 
of the state are fractionally charged (q = | for i^ = |) 
and obey fractional statistics (0 = ^ iox v ~ i). 

In Ref. 9, a similar scenario as the one outlined 
above for the FQH effect has been studied numerically 
in a ID system: These authors consider a topologically 
non-trivial flat band similar to t he mo del introduced by 
Su, Schricffer, and Heeger (SSH)i21iil at rational filling 
V — t: which is subjected to short-ranged interactions. 
Their numerical data indicates remarkable simila rities 
to the FQH setting. The quantized Berry phase^^E^ 



playing the role of the topological invariant of the non- 
interacting ID band-structure assumes fractional values. 
The system exhibits a ground state degeneracy of three 
when periodic boundary conditions are applied. The ele- 
mentary excitations of the system carry fractional charge. 

In this work, we present an exact solution of the 
model for the one-dimensional fractional topological 
phase (IDFTP) discussed in Ref. [9] in the physically 
relevant limit of a large band-gap where a projection 
onto the partially filled lower band is justified. This 
approach is analogous to the widely used projection 
onto the lowest Landau level in the FQH case. Our 
solution affords an intuitive physical interpretation of 
all the mentioned peculiarities of the IDFTP and allows 
us to scrutinize the key differences between the IDFTP 
and the FQH scenario at an analytical level. Even the 
phase diagram resulting from the competition between 
a nearest neighbor (NN) and next-nearest neighbor 
(NNN) interaction can be precisely understood at the 
level of our exact solution. 

In agreement with the general relation between 
Berry phase and entanglement entropy!^, we find that 
a fractional von Neumann entropy characterizes the 
reduced density matrix of a translation-invariant ground 
state of a bipartite IDFTP. However, when calculating 
the Berry phase of the interacting ground states we find 
a remarkable difference to the FQH case: In Ref. [HI it 
has been shown that the Hall conductance of a gapped 
system is insensitive to twisted boundary conditions 
(TBC). More precisely, the Hall conductance can be 
expressed as a constant Berry curvature defined on the 
torus of twisting angles. For the v = I: FQH case, 
Niu et a\W^ have shown explicitly, that the fraction- 
alized Hall conductance can be viewed as the Chern 



numbeiSEl over the enlarged torus of twisting angles that 
encompasses all three degenerate ground states. In the 
ID system under investigation in this work, the Berry 
phase represents the charge polarization of the systenP^l 
and is defined in terms of the Berry connection rather 
than the curvature which brings about a certain basis 
dependence. More precisely, we can find a basis of the 
ground state manifold in which only one state carries 
the total Berry phase of tt and the other two states are 
independent of the boundary conditions. However, in 
a translation-invariant basis, each ground state can be 
assigned a Berry phase of ^. This observation is closely 
related to the well known fact that a charge polarization 
only has a relative meaning depending on the choice of a 
reference unit celP^l. In contrast, the Hall conductance 
of a 2D system is a directly observable quantity with an 
unambiguously defined value. 

Furthermore, we investigate the entanglement spec- 
trum of the IDFTP employing the so called particle 
cu1p3 and find remarkable differences to the Laughlin 
z/ = I state^l in the thin torus limilpSl. In the quan- 
tum Hall case, the number of 'entanglement levels' (or 
equivalently, the rank of the reduced density matrix), is 
given by the number of ground states of the Hamiltonian 
for the Laughlin state itself, with a reduced number of 
particles, but at the original number of flux quanta^^. 
Thus, the particle entanglement spectrum of the ground 
state probes the quasi-hole excitations. Using the ex- 
act solution, we show that for the IDFTP the rank of 
the reduced density matrix is in fact much smaller than 
the number of quasi-hole states, and in this case merely 
probes the fermionic nature of the particles. 

Another important difference to the 2D case is that 
topological order in ID is always symmetry protected as 
long as particle number is conserved, see for instance, 
Refs. [2DH111 We explicitly show the symmetry protected 
nature of the IDFTP by constructing a gapped inter- 
polation to a trivial charge density wave without any 
polarization, involving the breaking of the protecting 
chiral symmetry. During this interpolation the total 
Berry phase of the ground states changes adiabatically 
from TT to zero. In contrast, an interpolation preserving 
the protecting symmetry involves a phase transition at 
which the Berry phase jumps to zero. 

We note that the experimental study of the predicted 
phase diagram for the IDFTP should be feasible in a 
synthetic system of cold atoms in an optical lattice. As 
has been demonstrated very recently, such settings even 
afford an experimental access to Berry phases^S 



SSH in 197^12111]. The tight binding Hamiltonian on 
which the IDFTP is constructed reads (see also Ref. [9]) 



Ho = ^ Cjdj+i + dLi c 



'j + 1^3' 



(1) 



where we have chosen unit lattice constant and Cj , dj 



are 



the annihilation operators of the two orbitals at site j. 
The Bloch Hamiltonian of this model can be written as 



/lo(fc) = v'-a-i 

v^ — cos(fc), v^ — sin(A:), u^ 



0, 



(2) 



where the ai denote Pauli matrices in the band pseudo- 
spin space. The spectrum of this Hamiltonian can be con- 
veniently obtained by taking the square E'^ = |v(fc)p = 
1, i.e., the model has completely flat bands. Furthermore, 
the Hamiltonian anti-commutes with the chiral symme- 
try operation cr^. The chiral symmetry can be viewed 
as the combination of a particle hole symmetry (PHS) 
operation C ~ a^K and the pseudo time reversal sym- 
metry (TRS) operation T = K, where K denotes com- 
plex conjugation. If one of the bands is filled, the system 
is gapped and its topological invariant S, is given by the 
winding number of the map k i— >■ (w^(fc), w^(fc))^ around 
the origin of the xy-plane^*. For our model which is a 
tight binding analog of the SSH modeP^I^, we obtain 
^ = 1. The physical consequence of this topol ogical ly 
non-trivial structure is the quantized Berry phasi 



ogicai t 



Lp = — ^ 



2t: 



dfc A{k) = TT (mod27r), 



(3) 



where A{k) — {lk\dk\lk) is the Berry connection of the 
Bloch states \lk) of the lower band. The Berry phase 
is related to the charge polarization P by 27rP = ip^ , 
see Ref. [12l The topological quantization is symmetry 
protected, i.e., if we allow for v^ ^ vfe can adiabati- 
cally connect phases with zero and tt Berry phase. The 
quantized polarization of the SSH model is illustrated in 
Fig. fll The noninteracting model in Eq. (pi) is readily 





FIG. 1. Trivial unpolarized model with v^ = 



= / 



v'" (top) with orbitals localized on the sites. Non-trivial SSH 
model with polarized orbitals localized on the bonds between 
two sites (bottom). 



II. MODEL AND ANALYTICAL SOLUTION 

We consider a two band lattice model similar to the 
dimer model for polyacetylene originally introduced by 



diagonalized as 



±, 



(4) 



where 7j,± = -75 (cjidj+i) are the annihilation operators 
of the eigenstates. From now on we consider a fraction- 
ally filled lower band and take into account two interac- 
tion terms, Vi for NN and V2 for NNN interactional 



Hj = Vi 



<i,j> 



niUj + V2 



<:i,j:s> 



(5) 



where rij — c^Cj + djdj. 



We are interested in the low energy physics coming 
from the interplay of the macroscopic degeneracy of 
states at the Fermi energy Ep = — 1 and the interactions 
described by Hj. Hence, we consider the non-interacting 
band gap of 2|w| as infinitely large compared to the scale 
of the interaction energies Vi, V2, i.e., a lowest band pro- 
jection (LBP). This assumption is similar to a projection 
to the lowest Landau level familiar in FQH physics. From 
a topological point of view there is a key difference be- 
tween continuum models (e.g. Landau levels of a homoge- 
neous electron gas in a perpendicular magnetic field) and 
lattice models. In the Landau level problem, the same 
Hall conductance of one quantum of conductance can be 
assigned to all Landau levels. Therefore, the topologi- 
cal defects of the Landau levels only add up to a larger 
and larger Hall conductance if more and more of them 
are filled. The physical reason for this is that all eigen- 
states of the underlying Hamiltonian describe cyclotron 
motions with the same chirality that is fixed by the di- 
rection of the magnetic field. In any lattice model, the 
situation is fundamentally different. The topology of a 
lattice model is defined relative to the total Hilbert space 
spanned by all bands. In more technical terms, the total 
Hilbert space of all bands can be seen as an embedding 
space for the bundle of occupied states. Hence, as far as 
the topological invariant associated with a given energy 
gap in a lattice Hamiltonian is concerned, the union of 
all bands above the gap is always the 'topological com- 
plement' of all states below the gap. Along these lines, 
we argue that all topological features of our lattice model 
must be encompassed by the LBP approximation. This is 
because the total Hilbert space of the two bands is topo- 
logically trivial so that mixing of the lower band and its 
'topological complement', i.e., the upper band can only 
perturb a topological state found in the LBP rather than 
leading to a richer topological structure. 

Within the LBP, we rewrite Hj in terms of the 
7i,± operators and discard all terms that contain 7+- 
contributions representing excitations of the upper band. 
Up to constant energy shifts, the interaction terms then 
read 



j 
^ nj{hj+i + 2nj+2 + fij+s) , 



4 

V2 



(6) 



where n. 



J 



At filling V = ^ a. state with the occupation num- 
ber pattern 001 obviously annihilates the Vi-term which 
immediately gives the gapped ground state at V2 = 

and explains its threefold degeneracy on the ring which 
has been numerically observed in Ref. |9l The de- 
generate states are obtained by translating the pattern 

001 (twice) by one lattice site and thus have the pat- 
terns 001,010,100. The size of the gap is ^ as follows 
immediately from the projected form ([6]) of the interac- 
tion Hamiltonian. This state is illustrated in Fig. [2j At 
filling ly = J, the state 0001 annihilates both terms in 
Hi and is ground state with a fourfold degeneracy on a 
ring and a gap ■^. 




oj — 7: 7- _ are the lower band occupation num- 
ber operators. 



FIG. 2. Gapped ground state at filling i^ ~ ^ and Vi > V2. 
Every third bond is filled with a particle (green dot). 

A non-trivial situation arises at filling i' — ^ ii both 
Vi and V2 are non-zero. In this case, placing two defects 
of the form 0001 and one of the form 0011 at constant 
filling fraction into the ground state 001 might become 
energetically profitable since the 0001 string annihilates 
the V2-term whereas the 001 pattern does not. Simple 
counting of interaction energies tells us that this defect 
pattern saves ^ and costs ^ due to the 0011 part. 
Hence, putting such defects becomes favorable at a criti- 
cal strength V2 = Vi at which the gap closes and a phase 
transition occurs. We confirmed the appearance of the 
phase transition at exactly this point in the limit of a 
large band-gap numerically. 



III. QUANTIZED BERRY PHASE AND 
SYMMETRY PROTECTION 



As already mentioned above, the completely filled 
lowest band of the non-interacting model is characterized 
by a quantized Berry phase of tt. Emplo ying the method 
of twisted boundary conditions (TBCJ^^EH^ -^q would 
now like to analytically calculate the Berry phase of 
the gapped interacting phase obtained within the LBP. 
The notion of TBC can be intuitively understood in 
the following way. One considers the (arbitrarily large 
but finite) physical system under investigation as one 
unit cell of a fictitious super-lattice. The lattice sites of 
the original lattice are now internal degrees of freedom 
of the super lattice, i.e., orbitals constituting one 
super-site. Upon Fourier transforming the super-lattice, 
each of these orbitals picks up a constant phase e^'^^ , 
where X labels the super-cell and (j) is the super-lattice 
momentum. In a particle number conserving Hamil- 
tonian the phase factors of creation and annihilation 



operators cancel out except for hopping terms crossing 
the boundary of a super-cell. 

Let us first apply this program to the non-interacting 
Hamiltonian (fTl) with a super cell of L sites. The model 
only contains NN hopping. Hence, the only term switch- 
ing the super cell is the hopping between orbital L of cell 
X and orbital 1 of cell X + 1. The Bloch Hamiltonian of 
the super-lattice associated with Eq. (IT]) is hence given 

by 



' L-l 

Hosi^) = I Y. 4^J+i + ^-c- I + (e'^c[di -I- h.c. 



(7) 



This model is still readily analytically diagonalized by the 
operators jj±{ip) = -7^{cj ± dj+i), j = 1, . . . ,L — 1 and 

7l,±(0) = 71 (ci ±e"^di). Remarkably, only jl^± de- 
pends on the momentum variable <j) of the super-lattice. 
For pedagogical reasons we would like to calculate the 
Berry phase of this model now in two equivalent ways. 
First, we interpret Eq. (l7| as a Bloch Hamiltonian with 
the occupied bands \uj-{(j))) — 7]_(0)|O). The Berry 
phase is then readily calculated as 



93-^ =-i 



277 

tl d0(O|7L,-(0)907l-(0)|O)=7r (8) 



An equivalent way to do this calculation is to interpret 



the Slater determinant |*o(0)> = UjTj.-Wl^) as the 
many body ground state of the Hamiltonian (I? I at half 
filling. In this many body language the Berry phase can 
be written as 



R 



2-K 



if- =-11 d,^ (vi/o(<^)|a^|*o(0)> = 

/o 

p2lT 

-i d0 (O|7L,-(0)907l -(0)|O) = vr. (9) 

Jo 

The latter approach is more useful for the generalization 
to the interacting model. 



A. The interacting model 

We now calculate the Berry phase of the interacting 
model with the Hamiltonian H = H^ + Hj (see Eq. (fTl) 
and Eq. (Isl)). As shown in Section Ih] the three degener- 
ate exact ground states of the interacting model for V2 < 



1,2,3. 



Fi at 1/3 filling are \^i) = UflV 7],+, _ |0), I ^ 
We now again apply TBC. The interaction Hamiltonian 
Hj, in particular in its projected form (see Eq. (|6|), 
does not depend on the twisting angle (f) since it does not 



contain any hopping terms. Therefore, even in the pres- 
ence of TBC, the interacting model can still be solved 
exactly by just replacing 7j._ with 7j,-((/!') in the defi- 
nition of the ground states |^;). The resulting ground 
states I'I'i ((/))) are obviously 27r-periodic in <j) so that the 
Berry phase defined in total analogy to Eq. (lo]) is a well 
defined geometric phase. Explicitly, we get 



/o 



(10) 



Two of the ground states thus have zero Berry phase 
whereas one of them has a Berry phase of tt. This result 
affords a simple physical interpretation. When imposing 
the TBC we go to a super-lattice and the Berry phase 
now describes the polarization of the model in the super 
cell. Only 1^3) contains an electron which is delocalized 
over two super-cells, namely the one created by j}^ _■ 
Hence j^s) has a polarization of i, or equivalently, a 
Berry phase of tt. In contrast, the two other ground 
state are unpolarized at the level of the super-lattice 
description and hence do not contribute to the total 
Berry phase. If we choose a different basis of ground 
states, we can distribute the Berry phase differently over 
the three basis states. For instance, upon combining 
the ground states IvP/) into momentum states, the latter 
contribute equally to the Berry phase. However, the 
total Berry phase (p^ — ^; (pf is always quantized to tt. 



B. Symmetry protection 

To demonstrate the symmetry protected nature of the 
present IDFTP, we now explicitly perform a gapped in- 
terpolation between this model and the trivial atomic 
insulator with the non-interacting Hamiltonian Hq — 



J. 



7t, 



^ ricj — d'dj. The Bloch Hamiltonian of this band 

structure is h{k) = ctz, so it obviously breaks the chiral 
symmetry. If we subject this model to the same inter- 
actions Hj as our original model, we again get three de- 
generate ground states at 1/3-filling of the lower band. 
These states are separated from the rest of the spectrum 
by a gap V2 ■ We will therefore assume NNN interactions 
to be present. Since Hq does not contain any hopping 
terms this model is completely insensitive to TBC and 
all three ground states \^i{(t))) = ni=o~ '^L+zl^) have 
a zero Berry phase as expected. We now consider the 
gapped interpolation 



H{\)=Hq{\)+Hi 
Ho{X) = Vl - XHo + VXHo 



(11) 



We first note that the Bloch Hamiltonian of the non- 
interacting model, Ho{X), reads 

hxik)^v'a, (12) 

^^ = VT^Xcos(k), yy = ^/^^^sin(fe), v' ^ V\ . 
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FIG. 3. The low-lying part of the spectrum of the interacting 
Hamiltonian -ff(A), with interaction parameters Vi = 1/10, 
V2 = 1/20 for system size L = 15 and periodic boundary con- 
ditions. We note that the ground state is threefold degenerate 
for each value of A. 
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FIG. 4. The low-lying part of the energy spectrum of the 
IDFTP with L — 15 sites and F = 5 fermions. The interac- 
tion parameters are Vi = 1/10 and V2 = (left panel), and 
Vi = 1/10 and V2 = 1/20 (right panel). 



Thus, this model also has flat bands, because i?| — 
\v{k)\'^ — 1, as was the case for Ho. For both A = 
and A = 1, the Hamiltonian has a three-fold degenerate 
ground state, separated from the excited states by a gap, 
because we assumed V2 > 0. The interaction Hj does 
not depend on the interpolation parameter A, and we 
find that H{X) has a threefold degenerate ground state 
and a gap for all < A < 1. 

In Fig. [Sj we show the low-energy part of the spectrum 
for a system with L — 15 sites at filling i^ = 1/3, with 
interaction parameters Vi = 1/10 and V2 — 1/20. We 
chose these parameters to be small in comparison to the 
non-interacting part of the Hamiltonian, so that the LBP 
is a good approximation. Indeed, the gap for A = 0, 
namely AE ss 0.02397 is close to the value Vi/A = 0.025 
valid in the LBP (see section lll| . For comparison, the 
gap for interaction parameters Vi = 1/10 and V2 = is 
AE w 0.02423. The gap for A = 1 is A^ = 0.05 = V2, 
the expected value. 

The plot confirms the existence of a gap throughout 
the interpolation. We only show the case of non-twisted 
boundary conditions = 0, but the spectrum is in fact 
gapped for all (j) throughout the interpolation. This is 
true, because twisting the boundary conditions only leads 
to a shift in the momentum k for Hamiltonians conserv- 
ing the number of fermions, and the one-particle energies 
do not depend on k. 

Finally, we note that we checked explicitly that upon 
increasing V2 from V2 = to V2 = 1/20 (while keeping 
Vi = 1/10 and A = fixed), one does not close the 
gap. The momentum resolved spectra for the two cases 
V2 = 0, 1/20 are displayed in Fig. |4) 

To calculate the total Berry phase, we changed the 
boundary conditions (j) from (j) — to 4> — 27r, in steps 
oi 6(j) = j^- We considered the system sizes L — 3 
(the non- interacting case), L = 6 and L = 9, as shown in 
Fig.[5]by the blue circles, red squares and green diamonds 
respectively. The interaction parameters are again Vi = 
1/10, V2 = 1/20. In the non-interacting case, the Berry 



phase can be calculated analytically. Namely, (p^ is given 
by half the solid angle mapped out by the curve of v{k) 
on the unit sphere when k changes from k — to k — 27r. 
This leads to ip^ — n{l—\/X), which is shown as the black 
line in Fig. [5J We find that the numerically obtained 
values using exact diagonalization for ip^ at i = 3 agree 
perfectly with the analytic result. The Berry phase in 
the interacting case for the larger system sizes with Vi = 
1/10, V2 = 1/20 deviates only slightly from the non- 
interacting result, confirming that the LBP is a good 
approximation in this regime. 

The Berry phase calculation clearly shows that the 
quantization of the Berry phase is protected by the chi- 
ral symmetry. Upon breaking this symmetry, we can 
smoothly change the Berry phase from ip^ — tt to 
if^ — 0, without closing the gap in the spectrum. 

It is also possible to interpolate to a trivial atomic 
insulator, without breaking the chiral symmetry. To 
this end, we introduce the non-interacting Hamiltonian 




FIG. 5. The numerically calculated Berry phase for the in- 
teracting Hamiltonian H{X) with Vi = 1/10, V2 = 1/20 as 
a function of the interpolation parameter A for L = 3 (blue 
circles), L = Q (red squares) and L = 9 (green diamonds). 
The black line is the analytic result for the non-interacting 
case, valid for L = 3. 
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FIG. 6. The map k <^ (v'' {k) , v^ (k)) of Eq. ([T4|, for various 
values of A: A = (black line), A = 1/5 (red dashed), A = 1/2 
(green dotted), A = 4/5 (blue dash-dotted), A = 1 (black dot). 
The black cross marks the origin, corresponding to iJ = 0, 
which lies on the curve with A — 1/2. 



HI, 



"^jCjdj + dljCj, such that the corresponding full 
interacting Hamiltonian reads 



H'{X) = Hl,{X) + Hj 



(13) 



il^(A) = VT^XHo 



rxH', 



The Bloch Hamiltonian of the non- interacting part Hq{X) 
now reads 

h'^{k)^v'a, (14) 

v'^ = VT^^cos{k) + \/A, yy = \/r^^sin(fc), w^ = . 

Indeed, the chiral symmetry is preserved. The spectrum 
of the Bloch Hamiltonian h'^{k) now depends on the in- 
terpolation parameter A and the momentum fc, namely 
El = 1 + 2,yX{l- X)cos{k). In particular, for A = 1/2 
and k — IT the gap closes due to a level crossing. We 
show the map k i— ?> {v^{k), u''(A;)) as defined by Eq. ( 14) 



for several values of A in Fig. [6j The Berry phase ~if^ 
corresponding to the lower band is given by tt times the 
winding number of the map k M- (v^ (k) , v^ (k)) around 
the origin. We find that for A < 1/2, the Berry phase is 
(fi^ = TT, while for A > 1/2, we have ip^ = 0. Precisely 
for A = 1/2, the winding number is not defined, signaled 
by the closing of the gap. 

For the interacting Hamiltonian H'{X), the scenario is 
very much the same, with the difference that the tran- 
sition occurs at different values of A. We determined 
numerically that the transition occurs at A w 0.34 for 
L = 6 and at A ~ 0.27 for L = 9. Due to its limited 
relevance for our main line of reasoning we did not study 
the precise depence of the gap closing on the system size 
in more detail here. 



IV. ENTANGLEMENT SPECTRA 

To gain further insight in the nature of the fractional 
topological phase in one-dimensional flatbands, we con- 



sider both the entanglement entropy and the entangle- 
ment spectrum. In the study of two-dimensional topo- 
logical phases, the entanglement entropy gives insight in 
the type of topological phase which is realized. Following 
the work of Kitaev and PreskilP^, and Levin and WerP^, 
we note that entanglement entropy S'ont associated with 
dividing the system into two (real space) regions A and 
B reads S'ont = ctL — 7 + 0(1/L), containing a (non- 
universal) contribution which scales as the length L of 
the boundary, as well as a universal constant 7 = log P. 
The quantity T) is the total quantum dimension associ- 
ated with the topological phase, which is a measure of 
the particle content of the topological phase (see, for in- 
stance, Ref. [28|) . In the context of the fractional quan- 
tum Hall states, this universal term has been determined 
numerically^^, and agrees reasonably with the theoretical 
predictions. 

We note that there are several interesting ways in 
which one can divide (or cut) t he sy stem into two parts 
A and B. The cut considered irpSHZlis the so-called 'real 
space' cut. In the 'orbital' cut, one divides the one- 
particle orbitals into two sets A and B. In the 'parti- 
cle' cut, finally, one divides the particles into two groups, 
with Na and Nb the number of particles in subsystem A 
and B, respectively. These different ways of dividing the 
system into two pieces probes different properties of the 
system, even though it might not be possible physically 
to actually perform the cut. 

The notion of the 'entanglement spectrum' was first 
considered in the context of the quantum Hall effect by 
Li and Haldane^. In short, the entanglement spectrum 
corresponds to the full spectrum of the reduced density 
matrix obtained by tracing out the degrees of freedom in 
part B. In comparison, the entanglement entropy com- 
bines all the eigenvalues into a single number. 

In this paper, we focus on the orbital and particle cuts. 
For these cuts, the entanglement spectra for fractional 
quantum Hall states were considered for various geome- 
tries in previous literatur^^*^'^^'^". 

Before we turn our focus on the IDFTP we study in 
this paper, we briefly mention some results concerning 
the particle entanglement spectrum. We focus on quan- 
tum Hall states for which a model Hamiltonian is known. 
This includes the Laughlin states, as well as many non- 
Abelian quantum Hall states, such as the Moore-Read 
statd^. The reason we focus on quantum Hall states 
with a known model Hamiltonian is that one can typi- 
cally obtain the number of zero energy ground states of 
these Hamiltonians, for an arbitrary number of electrons 
iVe, and an arbitrary number of flux quanta N^. When 
we divide the electrons into two groups A and B, and 
trace out the electrons in group B, we are left with a 
system in which the number of particles is reduced, but 
the number of flux quanta is unaltered. In the case of 
model Hamiltonians, one can show easily^ that the rank 
of the density matrix is bounded by the number of zero 
energy ground states of the model Hamiltonian, with N^ 
electrons but with the original number of flux quanta N^. 



It has been observed numerically that for the model 
quantum Hall states, this upper bound is indeed reached 
(see, for instance, Ref.[TO|l. Proving that the upper bound 
is reached has turned out to be hard and at presence a 
proof is only known for the Laughlin states^. 

In the case of the IDFTP, we perform a similar analysis 
by comparing the rank of the reduced density matrix to 
the number of ground states for a system with a reduced 
number of fermions, but with the same number of sites 
(playing the role of the number of flux quanta N^). 

The rank of the reduced density matrix, or equiva- 
lently, the number of levels in the particle entanglement 
spectrum, has been used beyond the realm of the frac- 
tional quantum Hall effect. In particular, it has been used 
to argue for the existence of so-called two-dimensional 
'fractional Chern insulators', see^for an early reference. 
It was found that the particle entanglement spectrum 
exhibits a 'gap'. The number of states below this 'entan- 
glement gap' was found to be given by the expected value 
from the quantum Hall states, showing the relation be- 
tween the fractional Chern insulators and the fractional 
quantum Hall states. 



A. The orbital cut 

We briefly discuss the entanglement entropy associated 
with cutting the system into two pieces. Because we con- 
sider periodic boundary conditions, we effectively cut the 
system in two locations when we trace out orbitals in sub- 
system B. 

To calculate the entanglement entropy, we start 
by recalling that if we work in the LBP, i.e., in 
terms of the fermions 7j,_ introduced in Sec. Hj the 
ground states are simple Slater determinants \Wi) = 

Hj-^o" T3j+i-\^)^ ' = lj2, 3. Thus, if we perform the 
cut in terms of the orbitals deflned by the fermions 
7j^_, the entanglement entropy 5*0111 = 0, because the 
ground states can be written as a single product |*() = 
|^j)a(8> \'(Pi)b- 

However, it is more relevant physically to consider the 
orbitals associated with the original fermions Cj and dj. 
In terms of these orbitals, the ground states are not sim- 
ple product states, and we will see that the entanglement 
entropy is nonzero. We first assume that only NN inter- 
actions are present, i.e. V2 = 0. Then, the form of the 
ground states in terms of the fermions Cj and dj is ob- 
tained from the explicit form of the operators 7j,- , which 
are given by 7j__ = ^(cj — ^j+i)- We divide the system 
into two subsystems A and B. If this division is such that 
none of the occupied 7j__ fermions 'has a component' in 
both system A and B, the entanglement entropy will still 
be zero, because we can still write |*/) — |V';)a (8) |V'i)s- 
If, however, the division is such that one jj^- fermion has 
a component in both A and B, the entanglement entropy 
will be given by ^cnt = In 2. Finally, if the cut is such 
that both boundaries contribute, we find ^ont = 2 In 2 
instead-^. 



We now briefly comment on what happens if we in- 
crease the interaction parameter V2, without making the 
assumption that we are working in the LBP. We will con- 
sider the case that the parameter V2 is small enough, 
such that we are still in the same phase as for V2 — 0. 
In this case, the ground state will have contributions in 
the upper band, and is more delocalized in comparison to 
the case V2 = 0. This means that upon cutting the sys- 
tem, the ground states will have longer range correlations 
across the boundaries, leading to an increased entangle- 
ment entropy, and an increase in the rank of the reduced 
density matrix. We confirmed this behavior numerically. 



B. The particle cut 

We pointed out above that in the case of quantum Hall 
states, the rank of the reduced density matrix, associated 
with dividing the electrons in two groups, is related to 
the number of ground states for the reduced number of 
particles, but with the same flux. 

We therefore now first take a quick look into the num- 
ber of ground states of the IDFTP, in case that the fill- 
ing is reduced from u = 1/3. To deduce the number of 
ground states of the IDFTP, with only nearest-neighbor 
interactions Vi present, we use exactly the same argu- 
ments as in Sec. [HJ We consider a system consisting of L 
sites, filled with F fermions. We assume that the filling 
1^ = J- < h, such that that there are ground states with 
the same energy as the ground states at filling v — 7;. 
The number of ground states is given by the number of 
ways in which we can distribute the F fermions over the 
L sites, such that no three consecutive orbitals contain 
more than one fermion. This number of ground states 
is precisely equal to the number of quantum Hall states 
on the torus, in the so-called 'thin-torus' limilP3n2Sl. Em- 
ploying some combinatorics shows that the number of 
ground states is given explicitly by p( p_i^), if F > 1 
(see also Ref. f37|l . For F = 0, there is trivially only 
one state. We checked numerically that this indeed gives 
the correct number of ground states of the IDFTP, when 
only NN interactions are present. 

Because we have the exact solution of the IDFTP 
model at our disposal, we can determine the particle en- 
tanglement spectrum exactly. As was the case for the 
orbital cut, it is easiest to do so in terms of the eigen- 
statcs I*,) = ntio"Slj+z._|0), I = 1,2,3. These states 
constitute a basis for the threefold degenerate ground 
state of the model at filling v — ^. Because we divide 
the system into two parts A and B, by means of dividing 
the fermions into two groups, we can perform the calcu- 
lation of the particle entanglement spectrum directly in 
the basis of 7j._ fermions. Rewriting these fermions in 
terms of the original Cj and dj fermions is, for the present 
cut, only a local transformation, which does not change 
the spectrum of the reduced density matrix. 

We concentrate on the ground state 1^3) = 



YijLi lli -|0)i which is the Slater determinant, such that 
each third orbital is filled. We now number the fermions 
and declare the fermions numbered 1, . . . , Na to belong 
to part A, while the remaining fermions belong to part 
B. Because the ground state j^a) is a single Slater de- 
terminant, it follows that the rank of the density ma- 
trix is given by the number of ways one can divide the 
Na fermions over the orbitals which are occupied in the 

original ground state 1^3) = rijiiTsi-lO) i'^^ refer to 
Refs. T^ and '55' for more details on calculating the par- 
ticle entanglement spectrum). 

It follows that the number of non-zero eigenvalues of 
the reduced density matrix is given hy [J ) = (^^ ) . 
Moreover, all these non-zero eigenvalues are equal to one 
another. We note that if we had calculated the reduced 
density matrix of the momentum eigenstates formed with 
I'^i), with I ~ 1,2,3, we would have found that the 
reduced density matrix has 3(jy ) degenerate non-zero 
eigenvalues. 

Having obtained the entanglement spectrum, we can 
now compare the number of entanglement levels to the 
number of ground states of the Hamiltonian at the re- 
duced number of fermions. Making use of the formula 
we gave above, we find that the number of ground states 
of the Hamiltonian with L — 3F sites, with Na > 
fermions, is given by I^C^^'"!!"') (or 1, if Na = 0), 
which should be compared to the rank of the reduced 
density matrix (^ ). We find that for Na = 0, these 
numbers are both equal to 1 trivially. For Na > 0, we 
find that the rank of the reduced density matrix is lower 
than the upper bound coming from the Hamiltonian. In 
particular, for Na = 1, the former is given by F, while 
the latter is 3F. For Na = F, the former is 1, while the 
latter is 3. In general, the ratio of the rank of the reduced 
density matrix and the upper bound is much smaller than 
one-third. 

The fact that the rank of the reduced density matrix of 
the particle entanglement spectrum is much smaller than 
the upper bound coming from the Hamiltonian marks a 
striking difference with the quantum Hall case, for which 
this upper bound is in fact satisfied. Indeed, the 'cor- 
relations' present in the fractional quantum Hall states, 
which are probed by the particle entanglement spectrum, 
are of a more non-trivial kind. In the IDFTP, they 
merely signal the fact that the ground states are Slater 
determinants of identical fermions. 

We close this section by making the following re- 
mark. In determining the particle entanglement spec- 
trum above, we did not make use of the fact that the 
interacting fermions are occupying a band with a non- 
trivial topology. In particular, we can give exactly the 
same arguments for the model in which we interpolated 
the bands to the trivial atomic insulator. In that case, we 
obtain exactly the same particle entanglement spectrum. 
Thus, the particle entanglement spectrum does not, in 
the present case, distinguish between the topological and 
trivial cases. This is not surprising, because for the mod- 
els we study, the particle entanglement spectrum probes 



the fermionic nature of the particles in the model. 



V. CONCLUDING DISCUSSION 

We investigated the IDFTP, which was first considered 
in Ref. [9l We solved the model in the limit where the 
interactions are small compared to the band gap, by pro- 
jecting the model onto the lowest band. This projection is 
analogous to considering the quantum Hall effect 'in the 
lowest Landau level'. The exact solution explains the ob- 
served threefold degenerate ground state and allows for 
the determination of the phase diagram. Although the 
IDFTP shares certain features of the fractional quantum 
Hall states, there are also crucial differences. 

The IDFTP exhibits a threefold degenerate ground 
state (when considering periodic boundary conditions), 
and fractionally charged excitations, just as is the case 
for the V = \ fractional quantum Hall state. We consid- 
ered the interacting model in a flatband with non-trivial 
topology (as in Ref. IH]), as well as in a trivial flatband. 
By choosing the appropriate interpolation (i.e., without 
breaking the chiral symmetry), we showed that one can 
adiabatically interpolate the IDFTP from the topological 
to the trivial case. That way we showed unambiguously 
that both the ground state degeneracy and the fraction- 
ally charged excitations should not be viewed as emerging 
because of the topological band structure, but rather as 
consequences of the CDW physics describing both cases. 
We also demonstrated that the Berry phase changes con- 
tinuously from LpB = TT for the topological Hatbands to 
LpB = for the trivial flatbands, in the interpolation 
mentioned above. 

Using the exact solution, we analytically calculated the 
Berry phase associated with the degenerate ground state 
in the case of the topological flatbands. This calculation 
revealed that upon changing the basis for the threefold 
degenerate ground state, one can change the relative con- 
tribution from each state to the total (quantized) Berry 
phase ipB — TT- This basis dependence can be understood 
by realizing that in one-dimensional systems, the Berry 
phase is a measure of the charge polarization^^, which 
depends on the choice of the unit cell. Thus, one can not 
entertain the objective notion of a fractionalized Berry 
phase. In contrast, in the two-dimensional quantum Hall 
effect, the fractionalized Hall conductance is a physical 
observable that is directly accesssible experimentally and 
hence represents objective physical reality. 

The basis dependence of the charge polarization al- 
luded to in the previous paragraph also plays a role in 
the entanglement entropy associated with dividing the 
fermionic orbitals into two sets. We showed that in the 
basis of 7j._ orbitals, appearing in the exact solution of 
the model, the ground states are single Slater determi- 
nants. Dividing these orbitals into two sets gives a van- 
ishing entanglement entropy. In terms of the original 
orbitals in which the model is phrased, the electrons are 
delocalized, giving a finite polarization. Using the CDW 
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basis for the ground state, this in turn gives rise to a fi- 
nite entanglement entropy pin 2, where p = 0, 1, 2 is the 
number of fermions delocalized over the cut. 

We also considered the entanglement spectrum associ- 
ated with dividing the fermions themselves into two sets. 
In the context of fractional quantum Hall states (and 
fractional Chern insulators), the particle entanglement 
spectrum of the ground state is directly related to the 
excitations of the system. This is signaled by the rank 
of the reduced density matrix, which equals the upper 
bound set by the Hamiltonian of the system itself. In 
calculating the entanglement spectrum for the IDFTP, 
we found that the rank of the reduced density matrix is 
in general much lower, and can be understood completely 
by considering the fermionic nature of the particles. This 
implies that the ground state of the IDFTP does not con- 
tain the correlations necessary to provide full knowledge 
of the excitations of the system, in contrast to the quan- 
tum Hall case. 

The exact solution revealed that in the presence of NN 
interactions only, there is a three-fold degenerate ground 
state at filling v = ^. Upon adding a repulsive NNN in- 



teraction, one finds that there is a (fourfold) degenerate 
ground state at filling v = ^. Thus, by merely changing 
the range of the two-body interaction, we can change the 
filling at which the ground state occurs from having an 
odd denominator to one having an even denominator. In 
the quantum Hall case, one can not perform such a sim- 
ply change in the filling fraction, because the fermionic 
nature of the electrons does not allow for a fermionic 
Laughlin state at even denominator filling. 

The Moore-Read quantum Hall state^ does have even 
denominator filling fraction. The model Hamiltonian 
which has the Moore-Read quantum Hall state as its 
ground state is a three-body interaction. In the flatband 
models we considered in this paper, one can also consider 
three-body interactions. In fact, it is straightforward to 
construct an electrostatic three-body interaction in terms 
of the 7j^_ fermions, which has the expected six-fold de- 
generate ground state at filling i. Although this model 
shares some properties with its quantum Hall cousin, it 
exhibits the same differences as the IDFTP in compari- 
son to the Laughlin states. 
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